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Nomenclature 

area [m 2 ] 

cross-section area of a conducting thermal strap [m 2 ] 
convective heat transfer area of vapor cooling [m 2 ] 
pipe diameter [m] 
diameter of capillary pores [m] 
specific heat of magnet [J kg ' 1 K' 1 ] 
specific heat at constant pressure of helium [J kg ' 1 K' 1 ] 


A 

A s 

Ay 

D 

d 


un 


P 

«T) 

f 

g(T) 


temperature dependent helium II conductivity function [cm 

hydrodynamic friction factor, non-dimensional 

fE [smkg - 1 K _1 ] 

(T - T b ) 2 * _L_ (T - T b ) 


1 + 


hL 


convection heat transfer coefficient [W m ' 2 K' 1 ] 
boiling heat transfer coefficient [W m ' 2 K' 1 ] 
constant heat transfer coefficient [W m ' 2 K' 1 ] 
Kapitza’s heat conduction coefficient [W m ' 2 K' 1 ] 

^quench 

J y(T) dT 

2 

^amb 

J y(T) dT 


KW 3 ] 
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k thermal conductivity [W m ' 1 K' 1 ] 

K non-dimension parameter, p c p v (f L )' 1 AT^ 3 , expressing the ratio of 

external and internal convective heat transport in He II 
L length [m] 

m mass [Kg] 

m vapor mass flowrate [kg s' 1 ] 

Nu Nusselt Number, hL/k 

P pressure [Pa] 

q heat flux [W m' 2 ] 

Q heat power [W] 

Qhx rate of heat transfer in vapor heat exchanger [W] 

Qpeak peak heat flux through thermal straps at T amb while cooling down 
Qstrap heat transporting rate through strap to bath [W] 

R thermal resistance, AT/Q [K W' 1 ] 

s entropy [J mole ' 1 K' 1 ] 

t time [s] 

t CD cooling down time [s] 

Atj cooling down period from quench [s] 

Atj cooling down period from ambient [s] 

T absolute temperature [K] 

T 0 thermal strap low end temperature [K] 

T b helium bath temperature [K] 
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T amb ambient temperature [K] 

T in vapor heat exchanger inlet temperature [K] 

T ou t vapor heat exchanger outlet temperature [K] 

T m magnet temperature [K] 

Tquench m agnet temperature at quench [K] 

AT temperature difference [K1 

ATfacii temperature drop at the cold end of thermal strap, across the boiling 

film pq 


AT 


K 


S 


temperature drop on Kapitza’s conductivity [K] 
velocity [m s' 1 ] 

specific volume of vapor [m 3 mole' 1 ] 
normal fluid velocity [m s' 1 ] 
volume [m 3 ] 

linear coordinate (along thermal strap [m] 


yCD 


c m (T) 


k(T) (T - T b ) 


v N ", 


[s m g- 1 K- 1 ] 


Z A/A, + Nu* \) 


Greek letters 

a thermal diffusivity, k/(pc) [m 2 s' 1 ] 

P geometric factor 
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5 thickness of a cylindric shell [m] 

o 

t \ viscosity [N s m' ] 

X, latent heat of evaporation [J kg' 1 ] 

p density [kg m' 3 ] 

ji chemical potential [J mole 1 ] 


Subscripts 

amb ambient 

av average 

fe fountain effect pump 

m magnet 

n normal fluid 

min minimum 

s super fluid 

st thermal strap 

svp saturation conditions 


v 


vapor 
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INTRODUCTION 

1.1 ASTROMAG - General Review 

ASTROMAG is a planned Particle Astrophysics Magnet Facility. Basically it is a 
large magnetic spectrometer outside the Earth’s atmosphere for an extended 
period of time in orbit on a space station. A Definition Team^ summarized its 
scientific objectives assumably related to fundamental questions of astrophysics, 
cosmology and elementary particle physics. 

Since magnetic induction of about 7 Tesla is desired, it is planned to be a 
superconducting magnet cooled to liquid He II temperatures. 

The general structure of ASTROMAG is based on: 

a. Two superconducting magnet coils. 

b. Dewar of liquid He II to provide cooling capability for the magnets. 

c. Instrumentation: magnet spectrometer 

1. Matter- Anti Matter Spectrometer (MAS) 

2. Cosmic Ray Isotope Spectrometer (CRIS) 

d. Interfaces to the shuttle and space station. 

Many configurations^ 5 '^ of the superconducting magnets and the dewar have been 
proposed and evaluated, since those are the heart of the ASTROMAG. Baseline of 
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the magnet configuration and cryostat as presented in phase A study (2) and the 
one kept in mind will doing the present study is shown in Figs. 1 and 2. 

Resupply is assumed to take place every 12 to 24 months assuming a cryostat 
design to hold sufficient helium for 2-5-4 years operation. 

ASTROMAG’s development schedule reflects the plan of launching to space station 
in 1995. 

1.2 Cooling System Specification 

Based on sources (2) and (3) the magnet coil cooling system must operate under 
five distinct conditions. These are in order of importance. 

1.2.1 Magnet charging 

1.2. 1.1 While the magnets are being charged, their temperature should be 
1.4 - 1.8 K. 

1.2. 1.2 Power dissipation in each of the coils will be approximately 250 mW. 

1.2. 1.3 Dissipation in the persistent switch will be approximately 2.0 W. 

1.2.2 Charged operation 

1.2.2. 1 The magnet coils should be at 1.4 -r 1.8 K. 

1.2.2.2 Parasitic heat input to each coil plus field decay totals approximately 


10 mW. 


Superconductln 



Fig. 1: Baseline ASTROMAG magnet cryostat. 1 
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COPPER BASELINE MAGNET PARAMETERS 


Number of Magnet Coils 2 

Number of S/C Layers per Coll 34 

Number of QB Layers per Coil 4 

Number of Turns per Layer 72 

Number of S/C turns per Coil 2448 

Number of QB turns per coil 288 

Coil Outside Diameter (m) 1,66 

Coil Inside Diameter (m) 1,56 

Space Between the Coils (m) 2.00 

Coll Width (mm) 184.00 

Magnet Seif Inductance (H) 33.52 

11 MJ Design Current (A) 810.09 

Coil Peak Induction (T)* 6.74 

Intercoil Tensile Force (kN)* 251# 

S/C Matrix Current Density (A/ sq mm)* 40 5 

Quench Energy at 1.8 K (micro-joules) 9.6 


* At the 11 MJ Design Coil Current 

# 25.6 metric tons 


Fig. 2: 


Copper baseline superconducting coil.^ 
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1. 2.2.3 Total stored energy is 11 MJ. 

1.2.2.4 Parasitics on the cryogen tank from other sources total approximately 
60 mW. 

1.2.3 Magnet quench 

1.2.3. 1 The magnetic field drops to zero in 2-5 seconds. 

1.2.3. 2 The coil temperature rises to 100-120 K. 

1.2. 3.3 Following a quench, enough He II should remain in the tank to allow a 
magnet cooldown and resumption of normal operation. 

1.2.3.4 The temperature in the cryogenic tank must not exceed the lambda point 
at any time while there is helium left in the tank. 

1.2. 3.5 The goal is a maximum temperature in the tank of 1.95 K. 

1.2.4 Cooldown from room temperature 

The coil cooling system should bring the temperature of the coils to below 2.0 K 
within 72 hours. 

1.2.5 Post quench cooldown 

The coil cooling system should return the temperature of the coils to below 2 K 
within 48 hours after the quench. 

1.2.6 As a design goal the dewar and the cooling system should be capable for 4 


years operation. 
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1.3 Scope of Study 

1.3.1 Thermal straps 

Detailed study of materials and sizes for steady state and transient analysis. 

1.3.2 Thermal straps and vapor cooling 

1.3.3 Hofmann’s Loop System 

That concept^ 14 ‘ 17 ^ is one of the ASTROMAG’s leading designs. Based on 
published analysis and experiments in addition to theoretical considerations, its 
main parameters to meet the requirements are derived. 

1.3.4 Suggestion: He II internal convection cooling 

A basic comparison is drawn between internal and external convection to evaluate 
the necessity of the last one. Demonstrating feasibility of the internal convection, 
a general description of the proposed system is derived. 
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2. STEADY STATE ANALYSIS 

2.1 Assumptions 

2.1.1 General layout and notation, on Fig. 3. 

2.1.2 There are two separate magnetic coils, each thermally coupled to cryogenic 
bath. 

2.1.3 The upper limit of bath temperature has to meet the required lowest 
temperature which is 1.4 K - 1.8 K at charging. 

2.1.4 Sharing temperatures drop, on Fig. 4. 
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Fig. 3: Steady state thermal model and notation. 




onKapitza’s temp, drop 
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Fig. 4: Sharing temperature drop on steady state. 
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2.2 Mathematical Model 
2.2.1 Fourier’s equation 


or 



Q = k dT 
A dx 


At steady state there is an almost linear profile because in the small range of AT = 
0.1 K the k is quite constant as shown in Fig. 5. 


dT _ AT 
dx ~ 


A _ Q 

L k * AT 


2.2.2 The highest heat load occurs during the charging mode. The dissipation in 
each coil is: 0.25 W. 

2.2.3 A higher value of A/L means that for a given cross section area, A, a 
shorter strap is required for the job, and v.v. 

2.2.4 The above A/L value is the minimal one (representing the longest strap). 
Higher values might be used. 
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2,3 Copper (RRR = 100) 

k(1.5 K) = 2.5 W cm' 1 K" 1 


A 

L 


0.25 W 


2.5 * 0.1 K 

cm K 


1.0 cm 


2.3.1 Example 

For a cylindric shell of D = 200 cm and thickness of 5 = 0.5 cm its length should 
not exceed: 

t 200 ^ 3.14 cm * 0.5 cm n. . 

L = = 314 cm 

1.0 cm 

For 5 = 2 cm, L = 78.5 cm. 

2.3.2 Example 

For one large coil instead of two with a double dissipation we get a maximum half 
length: 157 cm. 

2.3.3 Summary 

Under the steady state requirement it seems that there is no practical limitation 


on the minimal length of the thermal strap. 



23 


2.4 Aluminum (1100 F) 


k(1.5 K) = 0.15 


W 

cm K 


2.4.1 


A 

L 


0.25 W 


0.15 


W 


= 16.7 cm 


cm K 


* 0.1 K 


Example 

Cylindrical shell D = 200 cm, 5 = 0.5 cm, the maximum length 

L = 18.8 cm 


(It’s mass is 15.93 kg.) 

Stainless steel type 304 

k(1.5 K) = 0.0008 — 

cm K 


A 

L 


0.25 

0.0008 * 0.1 


= 3125 cm 


2.4.2 


24 


For a concept of supporting structural load by the cooling strap so that: 5 = 5 cm 
for a cylindrical shell of D = 200 cm. 


L _ 200 * 3.14 * 5 
3125 


= 1.0 cm 


(It’s mass is 25 kg.) 

2.4.3 Aluminum 5083 

k (1.5 K) = 0.015 _5L 
cm K 


A _ 0.25 

L 0.015 * 0.1 


1.67 cm 


2.4.4 ’High puritv aluminum RRR = 1500 


k(1.5 K) ~ 15 


W 

cm K 


_ - 0.17 cm 
L 
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2.5 Kapitza’s Thermal Resistance 


h k = 0.04 T 3 


W 

cm 


hk 



Q 

A * AT 


@ 1.5 K h k = 0.135 


For 


Q = 0.25 W and AT] 


0.25 W 


0.135 — — — * 0.1 K 
cm 2 K 


K 


Q 

h k * AT 


W 

cm 2 K 


= 0.1 K 


= 18.5 cm 2 


A cylindric shell of D = 200 cm has that cross section area for 5 = 0.03 cm. It 
means that a very small area is enough for providing a low Kapitza’s resistance. 
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2.6 Vapor Contact 

Under conditions of low gravity the He II might be distributed inside the dewar so 
that the cold end of the strap contacts vapor. In this case the heat transfer 
coefficient is of boiling. A low value of film boiling heat transfer coefficient^ is: 

h = 0.02 _W_ 

cm 2 K 


The worst case is when exactly the whole cross section area of the strap is exposed 
to vapor. If less is exposed then the h of boiling is replaced by that of Kapitza’s 
conductance which has a higher value. If a larger area is exposed to vapor then 
the equivalent heat transfer coefficient h A is bigger again. 


_ Q 

AT * h 


0.25 W 


0.1 K * 0.02 


W 

cm 2 K 


= 125 cm 2 
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2.7 Thermal Conductivity in Magnetic Field 

Thermal conductivity, more apparent for high purity metals, tends to decrease 
sharply under magnetic field. So far the available data are for A1 RRR = 1500 up 
to 4.0 Tesla and down to 7.0 K. (9) The present need for data goes down to 1.5 K. 

It is hard to extrapolate, however it is possible to draw the following trends. 

1. At low temp, (under 50 K) the k decreases with temp. 

2. For A1 RRR = 1500, at 7.0 K, at 4 Tesla, the k is reduced by a factor of 
8. 7/2.9 = 3.0, compared to 0 Tesla field. The change from 3 Tesla to 4 
Tesla at 7.0 K is only 3.02/2.93 = 1.03 (source No. 9 p. 152). 

3. The field dependence of k exhibits a saturational behavior which practically 
means that no more decrease occurs under higher field. 

4. Toward 0 K the k approaches exponentially close to zero values. Therefore, 
under topological reasoning of the B dependence of k we may conclude: 


for A1 RRR = 1500 at 1.5 K: 


k (4 Tesla) _ 1 
k (0 Tesla) 3 


At 1.5 K for H.P. A1 extrapolated from source No. 8 about Aluminum 1 is: 


k (0 Tesla, 1.5 K) = 15 


W 

cm K 


So assumably under previous conclusions: 
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k (7 Tesla, 1.5 K) = 5 


W 


cm K 


A | = 0.25 

L H.P.A1 5 * 01 


= 0.5 cm 


For H.P. annealed copper through extrapolation from 4 K to 1.5 K (source No. 8) 


k (0 Tesla, 1.5 K) = 30 


W 


cm K 


Assumably, again under the same degradation through magnetic field: 


k (7 Tesla, 1.5 K) = 10 


W 


cm K 


a. 


= 0.25 cm 


Annealed Cu 


Through sources (5) and (7) the ASTROMAG magnetic field is assumably about 7 


Tesla. 
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2.8 Mixed Case (Example) 

Assuming that the ends of the strap are to be of the structural alloy A1 5083 and 
its heart of H.P. A1 as shown on Fig. 6. 


By source ( ) for Al-Al electropolished contacts at 4.2 K > T > 1.8 K 

h 2 = 3.6 * 10 ' 3 T 2 ' 3 — 

cm 2 K 


Applying the above relation for 1.5 K 



h 2 


0.005 


W 

cm 2 K 


h l 


0.02 


W 

cm 2 K 


The altogether thermal resistance R 


R = 


Ki * A 


+ 


L 2 

K 2 * A 


l 3 

Ki * A 


1 

+ .I n 

hj * A 


2 

h 2 * A 


I W 


when 


high purity A1 
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Fig. 6: Example of a two material composed thermal strap. 
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R * Q = AT 


For our case: (Phase A design) 


R = 


15 40 15 A 1 

0.015 A + 5 A + 0.015 A + 0.02 

JL [1000 + 8 + 1000 + 50 + 400] = 


2 

A + 0.005 A 

2458 K 
A W 


R = 22 A = 2458 Q = 2458 * 0 25 = 6145 cm 2 
Q AT 0.1 

For D = 200 cm 8 = 9.7 cm. 

The main contribution comes from the high thermal resistance of the A1 5083 and 
the contact resistance. 

Conclusions 

1. The ends of the strap have to be as short as possible, like 3 cm each one. 

2. The contact thermal resistance has to be specially designed by enlarging the 


contact area. 
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Improved design: (Phase B) 

L x = L 3 = 3 cm A contact = 2 A 


R = — [200 + 8 + 200 + 50 + 200] = 


658 


A _ 658 * 0.25 

0.1 


1645 cm 2 


@ D = 200 cm, 5 = 2.6 cm 
Phase C 

Li = L 2 = 1 cm Aj = A 2 = 2 A A contact = 2 A 

R = _L [50 + 8 + 50 + 25 + 200] = 433/A 


A = 1083 cm, 


@ D = 200 cm 


5 = 1.7 cm 
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a. TRANSIENT ANALYSIS: 

COOLING DOWN THROUGH ONLY 


3.1 Assumptions 

3.1.1 General layout, on Fig. 7 

T(t) - temperature of the magnet as a lamped mass, and the high temperature of 
the strap. 

T 0 - the low temperature of the strap 

3 . 1.2 

T„ = T b + AT W1 
T b , bath temp. 

ATboi], temp, drop at the boiling film. 

3.1.3 Temperature profiles along the strap are qualitatively derived from the heat 
conduction factor (k) temperature dependence, as shown in Fig. 8. 

Assuming developed profiles or the quasi static case: 

q * f(x) 


k ilZ. = q = Const 
dx 
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Fig. 7: General layout and notation for transient analysis of magnet cooling 
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Fig. 8: Transient (qualitative) temperature profiles of magnet cooling down 
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dT = q_ 
dx k 


For low k, dT/dx is large and v.v. 

3.1.4 Developed temperature profiles 

The non-dimensional time parameter (at/1 2 ) ~ 0.5 can serve as a criteria for the 
quasi static (developed) case. 

a, thermal diffusivity, _!i- 

pc 


l 2 

t = _ 0.5 


a 


Assuming L = 15 cm, for copper RRR = 100 


@ 300 K: 


a = 


pc 


8.96 


0.42 


g 


cm 1 


W 

cm K 
* 0.38 


= 0.12 


cm 


gK 
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0.5 * 225 cm 2 
0.12 cm 2 /s 


= 937 s = 15.6 min 


For A1 (1100) 


@ 300 K: a = — = 0.86 cm 2 /s 

2.7 * 0.902 


t = 


0.5 * 225 

0.86 


s = 130 s ~ 2 min 


The cooling process is much longer (2.59 * 10 5 sec or 4320 min) so under that 
resolution of time it might be justified to concern it as a quasi static case. 

3.1.5 The straight line temperature profile serves in this case a guideline or a 
frame for the exact solution because: 

1. for k independent on T the developed solution is linear (for any value of 
constant k) 

2. over a wide range of T (100-300 K) the thermal conductivity is almost 


constant. 
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3.1.6 Boiling 


While cooling down the rate of heat removal is much higher than at steady state, 
the average power for 72 hours of cooling down is: 


26 * 10 6 J 
2.5 * 10 5 s 


100.4 W 


A typically low value of film boiling heat transfer coefficient is: 


h = 0.02 


W 

cm 2 K 


defined through: 


q = h * AT 


For the practical case of D = 200 cm and 5 = 0.5 cm 


9av 


100.4 W 
314 cm 2 


= 0.318 


W 


cm 


AT 


av “ 


q 

h 


0.318 


W 

cm 2 


0.02 


W 


15.9 K = AT^! 


cm 2 K 
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Assuming that the average heat flux is reached at 150 K, the AT av = 15.9 K is 
quite non-negligible. Notation shown in Fig. 9. 

Conclusion : 

The cooling down strap conduction should be solved coupled with the boiling 
process. 

3.1.7 The magnet might be concerned as made mainly of copper, though its mass 
of 650 kg and storing 11 MJ which raises its temp, to 120 K. 

For copper: 

120 

m f c dT = 650 10 3 g 17 £ = 11.05 MJ 
i S 


3.1.8 Total internal energy remove is derived through temperature extension to 
300 K 


300 


AE = m f c dT = 650 * 10 3 g 40 i = 26 MJ 

•i 0 


3.1.9 The mass of the strap is negligible compared to the mass of the magnet. 
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Fig. 9: Temperature drop across the boiling film and profile while magnet 
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3.2 Mathematical Model 

3.2.1 The heat flux at the hot end of the strap, q(L), represents the cooling rate of 
the magnet dT/dt. 


Pm c m 


= A q(L) = -k A \ — 
dt dx 


k=L 


The index m refers to magnet. 


3.2.2 Basic assumption of the model (discussed later) 


"dT "| , AT T - T 0 

Jc-l L L 


3.2.3 The ruling differential equation 


Pm c m 


(T - T 0 ) =0 
dt T v 01 


kA 


That is a first order, non-linear diff. equation. Its coefficients are temperature 
dependent. 


3.2.4 Through 3.1.2 and 3.1.6: 
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q(x=0) 

h 


T o = T b + AT boil = T b + 


Assuming quasi static (developed) case: 


q (x=0) = q(x=l) 


But: 


q(x=L) 


= k 


T - T 0 
L 


So: 


T 0 =T b 


. kfT-To) 

L h 


T 0 + — = T t 
0 Lh 1 


kT 

Lh 


T 


0 


T b 

1 


kT 

Lh 

k 

Lh 
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T - T 0 = T - 


T b - 


1 + 


kT 

Lh 

k 

Lh 


T -T b 

1 + Jl 

Lh 


1 _ 1 

fi ♦ Jl] 

1 

f x + J_1 

T - T 0 T - T b 

Lh j 

T - T b 

NuJ 


Nusselt Number 


Nu 


hL 

k 


^^boil - T 0 T b = 


Tk 

Lh 


- Tv 


Lh 


k 

Lh _ (T - T b ) 

i k 1 + N 
1 + 


AT boil = ( T " T b) 


Lh 
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3.2.5 From 3.2.3 


L 

A 


V m Pm c m 

fr - t 0 ) k 


dT 


-dt 


p m is a very weak function of T so: V m p m = Taking 1/(T - T 0 ) from 3.2.4: 


— m 
A 


c m (T) * 1 + 


k(T) 

Lh 


m 


(T - T b ) * k(T) 


J dT = - dt 


Let: 


Cm<D 


y(T) = 


i + m] 

Lh _ 


(T - T b ) * k(D 


— rn m f(T) dT = - dt 
A 


T 


quench 


/ 


2 


y(T) dT 



and corresponding Atj 


Let: 
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T 


amb 


J y(T) dT = I 2 


2 


and corresponding At 2 


Then: 


A _ m m * li 
L At x 


A _ m m * I 2 
L Ata 
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3.3 Minimal Cross Section Area for Heat Transfer 
We get: 


A amb 




(T - T b ) k 


dT = At 


The case of: 


k 

hL 


> 1 


represents one of the following situations: 

* L -> 0 

* k, very big 

* h, very small 

1 + JL — 

+ _ hL _ hL 


then: 


k 


k 


1 

hL 
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and: 


m m C c m 
Ah J T - Tb 


dT = At 


or: 


m 


m 


^min 


h * At 


T'amb 

; 


c m (T) 


T - Tv 


dT = 


m 


m 


h * At 


* Int 


Assuming that the magnet is represented by properties of copper 


300 


/ 


Cm(T) . dT = 0.55 i. 
T - 1.4 K 


Mmn 


6.5 * 10 5 g * 0.55 — 

5_ = 69 cm^ 


0.02 


W 


* 2.59 * 10 5 s 


cm 2 K 
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3.4 Model Solution 

3.4.1 General procedure 

a) Assuming L 

b) Through 3.2.5 solving (numerically) the integral of f(T) and getting A/L 

c) Getting A 

d) So, for every L we can find the A 

3.4.2 Copper (RRR = 100) 

From Table 1 and Fig. 10 we get the y(T) and Jy(T) dt. 


3.4.2. 1 . L = 5 cm 


II 


120 

f y(T) dT = 2.19 S cm 
{ S 


Atj = 48 h = 1.73 10 5 s 


300 

1 2 = / y(T) dT = 5.61 
2 


s cm 
g 


Ata = 72 h = 2.59 10 5 s 


(A) 

L 


v 


A 


m m ^1 
Ati 


6.5 * 10 5 g * 2.19. 
1.73 * 10 5 s 


s cm 


= 8.23 cm 


g 
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('A'! = = 6 5 » 105 » 5 61 cm = 14 cm 

J 2 At 2.59 * 10 5 


A = 



= 5 cm * 14 cm = 70 cm 2 


3.4.2.2 L = 30 cm (Table 1 and Fig. 10) 


Ii 


= 0.456 


s cm 
G 


h 


0.99 


s cm 
~G~ 



\ 

A 


6.5 * 10 3 * 0.456 
1.73 * 10 5 


cm = 1.7 cm 


(A) 

14 


6.5 * 10 3 * 0.99 
2.59 * 10 5 


cm = 2.48 cm 


3.4.2.3 L = 60 (Table 1 and Fig. 10) 


I 2 


0.54 


s cm 


g 


Table 1: Data and calculations for thermal straps transient analysis. 
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[ T . ( X ) 01 
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-1 


6.5 * 10 5 * 0.54 
2.59 * 10 5 


cm = 1.36 cm 


A = 60 cm * 1.36 cm = 81.3 cm 2 


3.4. 2.4 L = 100 cm (Table 1. Fig. 10) 


I 2 = 0.37 
2 g 


fA'l 

14 


= 0.93 cm 


A = 100 cm * 0.93 cm = 93 cm 


3.4.3 Aluminum 5083 

3.4.3. 1 L = 5 cm (Table 2. Fig. 11) 


h 


5.9 


g 
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(A) - 6 5 • 105 • 5 9 cm . 14.8 cm 
} 2.59 * 10 5 


A = L 


(A) 

14 


= 5 cm * 14.8 cm = 74 cm 2 


3.4.3.2 L = 20 cm (Table 2, Fig. 11) 

I 2 = 2.02 


(A-) 


6.5 10 5 2.02 
2.59 10 5 


cm = 5.07 cm 


A = 


7A^ 


•M 


= 20 cm * 5.07 cm = 101.4 cm 2 


3.4,3.3 L = 50 (Table 2. Fig. 11) 


h 


1.32 


s cm 


g 
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6.5 * 10 5 * 1.32 
2.59 * 10 5 


cm = 3.3 cm 


A = 


fA') 

14 


= 50 cm * 3.3 cm = 165 cm 2 


3.4.3.4 L = 80 cm (Table 2. Fig. 11) 


h 


1.15 


s cm 
G 



= 2.8 cm 


( A 


■\ 



= 80 cm * 2.8 cm = 224 cm 2 


A = L 


Table 2: Data and calculations for thermal straps transient analysis 
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10 4 y(T) 46.2 12.0 13.0 54.45 55.8 40.7 18.85 
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3.5 Constant Average Properties Solution 


L 

A 


-m 


1 

L_ 



k 


dT 

T -T b 


= -dt 


c m (T) c av 

k (T) -> k av 


L 

A 


“av 


“av 

hL 


w av 


dT 

T -T b 


= -dt 


L 

A 


' av 


f ^ 

i + 

hL 


a A X 


k av 


dT 

T -T b 


dt 

x 
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In (P - T„)|J = 1 

l 



T b 

T b 


_t 

x 


T b - 1.4 


T - 1.4 _ t 

0.6 t 
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3.6 Magnet Temperature Decrease as Function of Time 

Figure 12 shows the general behavior of the relative cooling down time as function 
of the temperature for different cases. 
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3.7 Starn’s Cold End Temperature 
3.7.1 General 
In 3.2.4 we have seen: 


T - T 0 = 


T - T> 


1 + 


hL 


So 


T 0 = T - 


T -T b 

i + Jl 

hL 



-T b 

TT“ 

hL 


3.7.2 Copper 
3. 7.2.1 L = 5 cm 


T (K) 

20 

40 

100 

160 

200 

260 

300 

^^9 

19.9 

39.7 

98 

156.2 

195.3 

254 

293 

EES 

0.1 

0.3 

2 

CO 

00 

4.7 

6 

7 
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3.7.2.2 L = 60 cm 


T(K) 

20 

40 

100 

160 

200 

260 

Bfl 

^^3 

18.1 

36.9 

80.5 

124.1 





1.9 


19.5 

35.9 

44.8 

58.1 



3.7.3 AI 5083 L = 50 cm 


T (K) 

10 

40 

80 

150 

300 

EES 

1.3 

9.4 

27.8 

65.9 

163 

T - T 

A A o 

8.7 

30.6 

52.2 

84.1 

137 


Figure 13 shows T 0 and T - T 0 for each T. 
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3.8 Heat Power 
3.8.1 General 

The highest heat flux occurs at the first instant of cooling down, T = 300 K 

q=-ki!Z Q=kAiI 

dx dx 

approximated through our model as: 

Q - k A Z—Zi =kA(T-TJ 

Li Li 

Qpeak = k * u (Iamb " ^o) 


3.8.2 Copper RRR = 100 
3.8. 2.1 L = 5 cm 

A = 14.0 cm T amb - T 0 = 7 K 

Li 


w 


* 14 cm * 7 K = 412 W 


Qpeak ~ 4.2 


cm K 
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From Chapter 3.1.6: Q av = 100.4 W 


Qpeak = 4 x 
Qav 


3.8. 2.2 L = 60 cm 

A = 1.36 cm T amb - T 0 = 67 

Li 


Qpe a k = 4.2 * 1.36 * 67 W = 382.7 W 


Qpeak 

^av 


3.8 


3.8.3 A1 5083 L = 50 cm 


A = 3.3 cm T amb - T 0 = 137 K 

L 


Oneak = 1.2 — 3.3 cm * 137 K = 542.5 W 

v cm K 


^ peak = 542.5 / 100.4 = 5.4 

Q av 
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3.9 Model Evaluation 

3.9.1 Exact integral solution 

The model ensures removal of the exact total internal energy. 

3.9.2 The model uses to some extent higher heat fluxes than the real ones. So to 
some extent it underestimates the real cooling down time. 

3.9.3 It better approximates the cooling range from 300 K to 120 K. In that 
range the thermal conductivity is almost temperature independent. So a 
longer portion of the temp, profile is linear. 

3.9.4 The ratio of maximal heat power to the average one is about 5, which seems 
acceptable. 
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4. TRANSIENT ANALYSIS: 

COOLING DOWN THROUGH CONDUCTION AND VAPOR 


4.1 Shortage of Sole Strap Conduction Cooling Down 

After quench the internal energy of 11 MJ has to be removed. That is equivalent 
to the magnet stored energy. The latent heat of evaporation for He II: 

X = 20 J/G 


So the required amount of He II for recovery would be: 


11 MJ 
20 J/g 


= 5.5 * 10 5 G = 550 kg 


equivalent to 3850 dcm 3 . That is about half of the intended tank’s volume. 




68 



Layout and notation of magnet cooling down by thermal strap and 
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4.2 Additional Cooling by Vapor 

By passing the boiled off vapor through the magnet more energy is subtracted 
from every unit of He II. Instead of utilizing only the latent heat before releasing 
the vapor it enables additional cooling by convection. Schematically it is described 
in Fig. 14. 

Vapor is generated according to the amount of heat transported through thermal 
strap to helium bath (dewar). That vapor cools the magnet while flowing through 
it. In the inlet vapor’s temperature is that of the helium bath and at the outlet, 
T out , according to the efficiency of the heat exchanger. 
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4.3 Short Heat Exchanger Case 
4.3.1 Mathematical model 


m m c m 


dT 

dt 


dT 


= -k^.A,-hA,(r-T b ) 


dx 


where A s , cross section of conducting strap, 

Ay, area of heat convection from magnet to vapor, 
T b , bath or vapor temp. 


Similarly to the transient model the dT/dx is approximated by (T - T 0 )/L 
m m c m ^ * Hi (T - To) + h A v (T - T b ) = 0 


T = 
1 o 


T - T> 


hL 


Assuming that the thermal conductance, h, is the same at the boiling film on the 
cold end of the strap and between magnet and vapor 
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m m c m — — + 


dT k A s T - T b 


dt L i + k 
hL 


h A V (T - 


dT kh A, „ _ . , . 

m m c m — rr = - ~ — — ~ " ^b) “ ^ 


dt hL + k 


m m c m 


dT h A g 


dt ^ + hL 


(T - T b ) - h A v 


— 



K 


_ B 


m 

Cm 

dT 


w 

i 

3 

1 

H 

QT^ 

As ... * A 

h 

— 


1 + hL 





k 



Tb) = 0 


(T - Tb) 


fT - T b ) 


It can be integrated similarly to the case of transient solution (without vapor) 
because the left wing is but temperature dependent. 
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4.3.2 Minimal convection and conduction areas 


m m c m 


dT = - dt 


h (T - T b ) 


hL 

k 


For L — > 0 or for k °° we get: 


m m c m dT 

IT (T - Tb) A s * A v 


= - dt 


(A 


s 



m m r c m 
h * At J T - T b 


m 


m 


h * At 


* Int 


That integral was calculated in 3.3.1, so that 


(A s + A v ) . = 69 cm 2 

V s v /min 


A s = A v = 34.5 cm 2 


For A w /A s = 1 
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4.3.3 Solution (example) 

Instead of integration we might have an approximate solution based on the 
previous transient case integration (without vapor cooling). Rearranging the diff. 
equ. to: 


-m 


1 + 


m 


hL 


m 


(T-Tb)k 1 + 


'l + 

V k 


dT = -dt 


J 


Let: Z = 


1 + 


1 + 


hL 


The present solution differs from the previous one by Z. 


Assumptions 

4.3.3.1 For copper RRR = 100 or A1 RRR = 1500 or other high conductance 
material, for L < 100 cm, h = 0.02 W/cm K 


hL 
— < 

k 


1 
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4.3.3.2 For A1 5083 and L = 20 cm on average along the whole temperature range 



4.3.3.3 


—I = Const 
A s 


4.3.4 Results 


1. For high k materials 



2. For A1 5083 and L = 20 cm 


Z = 


A v 

1 + 1.3 — 
A a 


3. The solving integral is proportional to the integral in the case of sole 


conduction. Denoting: 
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I 



fl + 

c m 

hL J 

(T 

-T b )k 


hL 

k 


dT = (II)i 


So: 



and 


m 


m 


m 


<id * - ^ 


At 


h 

Z 


4.3.4.1 Assuming A/Ag = 1 that: 

for case 4.3.3. 1 Z = 0.50 and 

for case 4.3. 3.2 Z = 0.43 

To sum up: 

4. 3.4. 2 For copper RRR = 100, A1 RRR = 1000 and k/Lh » 1 
L = 30 cm A s = 37.2 cm = A v 


L = 60 cm 


A s = 40.8 cm = Ay 
Ag = 46.5 cm = Ay 
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L = 100 cm 


4. 3.4. 3 For A1 5083 

L = 20 cm A s = 43.6 cm 2 = Ay 


4.3.5 Reduction in boil off while cooling down 

The parameter Z represents the ration between the heat flux in case of sole 
conduction and the case of conduction plus vapor cooling. 

Under the above assumptions, when Z is constant, it means that the fluxes are 
proportional at every moment to those of the case of cooling through thermal strap 
only. 


Therefore the total amount of energy removed through conduction in presence of 
vapor cooling is lower by the factor Z. 


The boil off is reduced by the factor Z. 
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4.4 Ideal Heat Exchanger 

The vapor passes through a heat exchanger in the magnet. It’s outlet temperature 
reaches the temperature T of the magnet. 

m m c m 4£ - -* c p (T - T b ) - (T - T 0 ) 

m, helium vapor mass flow rate 
c p , helium vapor heat capacity 

hi = Qstrap/^ 

Q s trap, heat power transferred through straps 
X, latent heat of evaporation 

Qstrap ’ ^ (T - T ») 


* k A T - T b 

XL, k 
1 + — 

hL 


Through chapter 4.3.1: 
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m m c m 


dT 

dt 


k A Cp (T - T b ) 2 - A 


X L 


1 + 


hL 


L , k 
1 + 

hL 


(T - T b ) 


L = m 


c m 


T A11 m 

A 

k c p 

/T - Thl 2 +- 

k 


f k ^ 

\ of 

r k ^ 

X 

1 + _ 


i + _ 


l 


hL J 


(T - Tfc) 


g(T) = 


■'m 


k c r 


hL 


(T - T„) s 


(T - T b ) 


hL 


for A1 5083 L = 20 cm according to Table 3 and Fig. 15. 


I 2 = 0.12 


s cm 

S 


A = 6.5 * 10 5 * 0.12 
L 2.59 * 10 5 


= 0.30 cm 


A = 20 cm * 0.15 cm = 3 cm' 


Table 3: Data and calculations for thermal strap and ideal heat exchanger transient analysis. 
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The peak heat flux through thermal strap: 

<W - k a Tamb r ' T ° ■ A fr amb - T b) 

L L i + J_ 


= 0.3 cm * 0.3 — (300 K - 1.4 K) = 26.9 W 

cm K 


Through chapter 3.1.6 we have seen: 

Q av = 100.4 W 
or 

Qpeak _ 26.9 _ n 9fi 
100.4 

- vfav 


We might wish to compare the peak heat flux through the thermal strap to the 
case of sole strap cooling (without vapor) for A1 5083 L = 50 cm which is 2.5 times 
longer than the present strap (with vapor cooling) we got in chapter 3.8.3. 

Qpeak _ 5 4 

Qav 
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Put another way, the peak heat flux is reduced by a factor of 5.4/0.26 = 20.7. It 
means that the perfect heat exchanger is making the major cooling work. 
However, the thermal straps are essential to transport heat for vapor generation. 
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4.5 Real heat exchanger 


m m * c m 


dT 

dt 


T - T 

= - kA_ tl 


- Q 


HX 


The rate of heat transfer in a real heat exchanger, Qhx, is ruled^ ^ by the 
logarithmic average of temperature differences: 


. t t * AT c - AT W 

Qhx = U * A v * — — 

In c 


AT 


W 


AT C *T - T b 

AT w = T - T oot 


Qhx = U 


* Ay 



-T b 
-T b 
“ T 0 ut 


(4.5.1) 


Energy balance on vapor stream 


Tout 




141 * c p fTout - T b) = Qhx = U ♦ Ay 
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T - T b _ U A v 
T - T out m * c p 


(4.5.2) 


f 


Vapor’s flowrate derives from heat power conducted to the bath through the 
thermal strap: 


m X = A 


* k 


T - T 0 
L 


Through chapter 4.3.1 


A * k T - T b 
X L ! k 
hL 


(4.5.3) 


The overall balance differential equation: 


rn m * c p *. 


dT 

dt 


_ (T - T b ) + U * A v Tput — = 0 


hL 


In 


T - Tv 


(4.5.4) 


T - T, 


out 


In its final form the model consists of: 
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(a) differential equation: 


dT 

dt 


kA 


r kl 

ITlm * C n * L 

1 + _ 

m p 

l hL J 


(T 



+ (T out - T b ) 


= 0 


(b) transcendental equation: 


In 


T - T b 

T ~ Tout 


U * X *L A v 
c p * k A 



_k_ 

hL 


j 


(T -T b ) 


Well approximated: 


U = h 


Through the transcendental equation (b) Tout is solved as function of T. Then 
equation (a) is integrated to get T(t). 
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5. THERMAL STRAP SUMMARY: 

CONCLUSIONS AND DESIGN 
REMARKS 

5.1 The constraints derived from steady state analysis and from transient 

analysis are plotted for A1 5083 on Fig. 16 and for Cu RRR = 100 on Fig. 17. 


5.2 The steady state constraints are much more restrictive, so they should 
dominate the design. 


5.3 For a steady state mixed case 


Ll L 2 1 

= + + 

k 1 A l k 2 A 2 h c A 3 


1 ^ AT 

+ < 

h b Ai Q 


R, thermal resistance 

W 

h c , contact heat transfer coefficient (A 3 ) 
h b , boiling heat transfer coefficient (A x ) 

5.4 As a general trend increasing \ enables reducing the A s for a desired after 
quench boil of amount. 


5.5 


Because of the domination of the steady state case, we may avoid the 
transient analysis for the mixed case (different material for strap). 
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5.6 The contact conductance between the metal seems to be significant so the 
contact areas are suggested to be enlarged. 

5.7 Using high purity metal (high conduction and high elongation) in between 
structural ends (low conduction and stiff) seems to be a mechanical 
advantage. It enables self adjustment, lower tolerances and easier assembly. 


pure conduction 
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Fig. 16: Restrictive sizing characteristics of aluminum thermal strap at various modes of operation. 
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[ z ] V 


Fig. 17: Restrictive sizing characteristics of copper thermal strap at various modes of operation. 
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6. PRELIMINARY DESIGN ATTEMPT 

6.1 Principles 

6.1.1 To meet the general layout and interfacing according to source 2 (Fig. 1, 
p. 3). 

6.1.2 In order to shorten the cold end of the strap (being made of low conductive, 
construction alloy) it will be part of the dewar design . 

6.1.3 Contact area enlargement 

6.1.4 Using A1 RRR = 1500 

6.1.5 Vapor heat transfer area if considered as the whole magnet surface. 

Ay = 166 cm * tc * 18.4 cm = 9590 cm 2 


6.2 Sizing 
D = 166 cm 
5 (contact) = 10 cm 
A1 5083 

L = 2 cm + 2 cm 


5 (conduction) = 2 cm 
H.P, A1 
L = 20 cm 


A s = 166 cm * 3.14 * 2 cm = 1043 cm 


A contact = 5220 cm2 
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6.3 Conduction 


R = 


4 

0.015 1043 


20 

5 1043 


2 

+ 0.005 5220 


1 

0.02 5220 


A1 5083 H.P. A1 contacts film boiling. 


R = 0.255 + 0.004 + 0.076 + 0.01 = 0.345 JL 

W 


AT = R Q = 0.25 W * 0.345 _ = 0.08 K < 0.1 K 

W 


6.4 Vapor Cooling 


= 1 
hb 


1 

0.02 * 9590 


0.052 


_K 

W 


R = 0.345 
R v 0.052 
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6.5 Thermal Stresses 

To ensure that the developed thermal stresses in strap between 300 K and 1.5 K 
will not cause failure. 

The elongation for H.P. A1 (to failure) is about 30%. The thermal expansion in our 
case is 

L293 ~ L ° - = 415 10' 5 = 0.415 % 

^273 

It seems to be able to withstand the tensile stresses. 

There is no problem of buckling since the stresses are tensile. 


Derived general layout shown in Fig. 18. 


He H DEWAR 
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Fig. 18: Preliminary design attempt of a thermal strap. 
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7. HE II INTERNAL CONVECTION COOLING SYSTEM 

7.1 Introduction 

Many space-based low temperature systems use stored He II (T < 2.2 K) to supply 
cooling. Included among these are infrared telescopes such as IRAS, COBE and 
SIRTF, as well as specialized systems like ASTROMAG, a particle physics 
experiment. In all cases, these cryogenic systems require a thermal link between 
the liquid helium dewar and the principal low temperature device, i.e., telescope 
or superconducting magnet. A number of possible approaches exist to make this 
thermal link. In the case of infrared telescopes, the main approach has involved 
direct thermal contact through support structure or high conductivity thermal 
straps. This method has been successfully applied to the IRAS and COBE 
cryogenic system. By contrast, the ASTROMAG cryogenic system has a leading 
design which involves the use of a self-activated fountain effect pumping system. 
The use of He II in direct contact with the superconducting magnets may be 
desirable because of the need to ensure a higher degree of thermal stabilization. 

In the present chapter we consider an alternative approach for cooling large He II 
cryogenic systems in zero-gravity. This approach involves the use of a kind of 
heat pipe which relies on counterflow heat transport in the He II contained in 
channels linking the dewar to the low temperature device. If it could be made to 
work effectively, a He II heat pipe system would provide better thermal contact 
than the high thermal conductivity straps. It also may be simpler to implement 
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and more reliable than the self-activating fountain effect pump. 

The self-activating fountain effect pumping system (Hofmann loop) has been the 

n 

preferred concept for cooling of the superconducting magnets in ASTRO MAG. 

The device uses a fountain pump activated by the parasitic heat leak into the 
magnet to force the He II to circulate through tubes surrounding the magnet . 14 
Fountain pumps have been developed for space-based systems as part of the 
Superfluid Helium On-Orbit Transfer (SHOOT) demonstration. Advantages 
include simplicity of design, where the pump consists of a heater and porous 
ceramic plug, and the tendency to suppress cavitation, a problem more common in 
centrifugal pumps with He II. The principal disadvantage, that of a low 
thermodynamic efficiency, is overcome in the Hofmann loop by relying on the heat 
leak into the system to drive the pump. 
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7.2 He II Internal Versus External Convection 

To consider whether the Hofmann loop makes sense for ASTROMAG, one must 
compare the relative merits, mostly improved cooling, with the increased 
complexity associated with the added components (pumps, heat exchangers and 
flowmeters) in the cryogenic system. Some of these issues require more detailed 
design; however, scaling arguments can be used for comparison. One issue to 
consider is whether the forced circulation will result in any benefit to heat 
exchange between the magnet and He II reservoir. Two factors are important 
here. First, the heat exchange process between the surface of the cooling channel 
and the He II is not affected by flow rate. This fact, which has been demonstrated 
experimentally 19,20 , results from the surface heat transfer being dominated by 
Kapitza conductance. Kapitza conductance is a phonon-mediated process which is 
not affected by the relatively small velocities of He II forced flow. 


Forced flow can have a significant impact on the heat transport in bulk He II. 
However, the question remains as to whether this process is of benefit to the 
ASTROMAG cooling system. One can determine the temperature profile in forced 
flow He n by analysis of the steady state He II energy equation 


1 dT^ m 

1 dx) j 


- p C V 


dT + d Vn = 0 


dx 


dx 


_d_ 

dx 


( 7 . 1 ) 
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where f = A p n /p s 3 s 4 T 3 , the He II thermal resistance parameter with A being the 
Gorter-Mellink mutual friction parameter. Solution to eq. (7.1) depends on the 
particular boundary conditions; however, it is possible to show by dimensionless 
scaling that the parameter 

K = p c p v (fL) 1/3 AT 2 ^ 3 (7 * 2) 

provides a measure of the degree to which heat is carried by forced convection 
compared to internal convection. If K' > 1, there will be a significant increase over 
static internal convection. Two variables in eq. (2) which must be determined are 
the fluid velocity and temperature rise AT in the helium. 

The He II flow rate through the system is determined by the performance of the 
fountain effect pump. For an ideal pump working against negligible pressure 
head the mass flow is given by 

m = Q (7.3) 

Sin Tin 


For a 10 mW heat generation at 1.4 K, eq. (7.3) yields a corresponding mass flow 
rate of about 50 mg/s. This mass flow rate results in a pressure drop through the 
hydraulic circuit. In turbulent flow, it has been shown that pressure drop in He II 
can be analyzed in terms of classical friction factor correlations. 21 For the sake of 
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simplicity, we assume that the Blasius correlation is appropriate for describing the 
friction factor and the pressure drop can be given by 


Ap 




(7.4) 


where f D = 0. 079/Re L/4 , which is appropriate for Reynolds numbers in the range 
10 4 . Substituting for the friction factor and rewriting in terms of mass flow rate, 
we obtain 


Ap 


0.241 __L_ 
D 19/4 



fQ 


V/4 



(7.5) 


The magnitude of the temperature change AT can be estimated by consideration of 
the fountain pump. If we assume that the pump has ideal behavior, then 
London’s equation applies, Ap = ps AT. The corresponding temperature difference 
is that which must be maintained across the pump to produce the required mass 
flow rate. This may not be the highest temperature in the flow loop, but it would 
be expected to be close except possibly at very high mass flow rates. Combining 
London’s equation with eq. (7.5) yields, 


AT = 0.241 


D 


19/4 


1/4 

% 

p 2 S 


Q_V/4 

sT 


(7.6) 


v j 


Ct5L 
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As an example, for a 10 mW heat generation at 1.4 K with L = 30 m and D = 7 
mm, eq. (7.7) predicts a temperature rise across the circuit of about 50 pK, while 
250 mW would result in AT ~ 14 mK. 


To determine the relative importance of forced convection heat transfer, we can 
compute the value of K' by substituting eq. (7.6) and (7.3) into (7.2). The result is 


K = 0.493 


q 13/6 C p L t ^ 6 

D 31/6 p4/3 g 17/6 T 13/6 


(7.7) 


As an example, for a 10 mW heat leak at 1.4 K and for ASTROMAG cooling loop 
geometry, eq. (7.7) predicts a K' ~ 5 x 10' 4 . During charge when Q = 250 mW, we 
obtain K' ~ 0.5. Note also that if the temperature increases, the value of K' 
further decreases because of the strong temperature dependence to the entropy. 
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7.3 He II Heat Pipe System 

The above result that forced convection contributes little to the heat transport in 
an ASTROMAG cooling loop except under extraordinary conditions calls into 
question the use of a self-activated fountain pump system in this application. The 
enhancement of cooling is small and the added complexity of design significant. 
Therefore, it seems appropriate to consider alternative cooling schemes for this 
system. One option is to use a simple static cooling loop which relies on internal 
convection to remove the steady state heat generation. We refer to this device as 
a He II Heat Pipe. Such a cooling system could be equipped with a fountain pump 
which would operate only for emergency conditions such as recovery from quench. 

A schematic representation of a He II Heat Pipe system for ASTROMAG is shown 
in Fig. 19. The magnet is indirectly cooled by a loop containing static He II. Each 
end of the loop penetrates the main helium dewar. Porous plugs located at the 
loop ends are needed to ensure that the helium within the cooling loop remains 
above saturation pressure. The dewar must also contain fluid acquisition devices 
(FAD) so that the porous plugs will be in continuous contact with He II. The loop 
is equipped with a heater (H) and bypass valve (BV) to allow the porous plug to 
function as a fountain pump under extraordinary operation such as magnet 
cooldown or recovery from a quench. 


There are two primary issues which must be considered when designing a He II 
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Fig. 19: Schematic representation of He II Heat Pipe system for 
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Heat Pipe system. These are that the loop must have sufficient cooling capability 
to remove the expected heat leak and that the helium within the loop must not 
cavitate. The heat removal capability of a cooling loop is determined by its 
physical dimensions and the allowable temperature difference within the system. 
This problem can be analyzed by knowing the total heat load on the loop and 
applying the heat transport equations appropriate for He II. 
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7.4 Temperature Drops Along the Heat Pipe 

For the purpose of analysis, we assume that the cooling loop consists of a one- 
dimensional channel of length L and diameter D, which receives a heat load Q 
distributed uniformly along its length. The loop is terminated at both ends by 
porous plugs, where we model as large number, N, of parallel capillaries of length 
1 and diameter d. The channel diameter and heat flux within the helium is taken 
to be large enough to obey turbulent heat transport so that the temperature 
gradient is given by 

£1 = f q 3 (7.8) 

dx 


where q is the heat flux W/cm 2 . For a channel of length L with constant 
temperature boundary conditions, the temperature difference between the middle 
and ends can be obtained by integration of eq. (7.8) to yield, 


AT X 


Q 


\3 



J 


f L 


(7.9) 


For typical ASTROMAG parameters, Q = 10 mW, D = 7 mm, L = 30 m and for T = 
1.8 K, f = 10' 13 m 5 K/W 3 , then ATj^ ~ 1 pK. Even under extraordinary conditions 
with Q = 250 mW the temperature difference is 13 mK, a generally small value. 
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The temperature difference across the porous plug is more complex owing to the 
normal fluid drag which contributes to the loss mechanism. Assuming that the 
porous plug can be modeled as N parallel capillaries, the temperature gradient is 
given by 


dT = M n q j. f q3 
dx (psd) 2 T 


( 7 . 10 ) 


The first term on RHS is the normal fluid drag. The normal fluid viscosity is % 
and P is a geometrical factor which depends on channel geometry but is 32 for 
circular cross-section. For a porous plug of a given porosity, a = N7td 2 /4Ap, and 
length 1, the temperature difference across it, AT 2 , may be calculated by 
integrating eq. (7.10). 


AT o = 


_g_ 

d 2 


f \ 

Q 

° A P j 


f 

Q 

v oA, 


1 


( 7 . 11 ) 


where we define g s p Tj n /(ps) 2 T, the laminar flow thermal resistivity function. 

Note that g - T‘ 12 so its contribution to the overall temperature difference depends 
strongly on the operating temperature range. Under most conditions of operation 
for ASTROMAG, the first term in eq. (7.11) will dominate the temperature 
difference across the plug. However, depending on the choice of pore diameter, 
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the turbulent term may have a significant effect. 

An approximate numerical estimate for the contribution that the porous plug 
makes to the overall heat removal of the He II heat pipe can be obtained by 
inserting typical values in eq. (7.11). At 1.8 K, g ~ 4 x 10' 15 Km 3 /W. For a steady 
heat leak of 10 mW, d = 10 pm, a = 0.5, 1 = 1 cm and Ap = 0.4 cm 2 , we obtain AT 2 
= 0.18 mK. Even for extraordinary conditions with Q = 250 mW, AT 2 ~ 5 mK, an 
acceptable value. A summary of different operating conditions and pore sizes are 


listed in Table 4. 
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7.5 Porous Media Evaluation 

The main purpose of the porous plug is to provide a pressure head so that the 
fluid within the channel will not cavitate when it receives a modest temperature 
rise due to heat exchange. In laminar flow, the pressure gradient is the result of 
normal fluid drag and can be written 


dp = P ^In v n 
dx d 2 


( 7 . 12 ) 


where v n is the normal fluid velocity. Equation (7.12) is simply the Poiseuille 
expression for the normal fluid component. In counterflow, v n = q/psT. Further, 
the pressure and temperature gradients are simply related through London’s 
equation. If we model the porous plug as a number of channels of length 1 
diameter d undergoing laminar flow, then the above equations can be integrated 

W 

to establish the size of the pressure difference across the plug for a given heat 
flux. 

Above the critical velocity, London’s equation is no longer valid owing to the onset 
of turbulence in the fluid. The effect of turbulence is largest in the temperature 
gradient because of the existence of the mutual friction term, second term in eq. 
(7.10). In some treatments of pressure gradient in counterflow He II, there is a 
second non-linear term which adds to eq. (12) such that 
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dp = PVn + o 00494 [Refrjl 1 - 75 (7.13) 

dx d 2 

which is appropriate for Re > 1200. Since the contribution by the second term is 
only significant for relatively high velocities, we will neglect it here and treat the 
temperature gradient explicitly in terms of eq. (7.12). This is the so called Allan- 
Reekie Rule. 

The situation we encounter with respect to the porous plug operation in reduced 
gravity is best described in terms of the liquid helium phase diagram, Fig. 20. We 
assume that the helium within the storage dewar is at saturation, thus being 
represented by a point b with a bath temperature T 2 . The helium within the heat 
pipe will be at a higher temperature, T lf as a result of heat absorbed from the 
residual heat leaks or charging losses to the magnet. If the porous plug behaved 
as an ideal fountain pump then the pressure difference would be given by 
London’s equation, which corresponds to a process of constant chemical potential, 
p = const. The pressure in the loop would be given by point a', which is well 
above saturation and therefore would prevent cavitation. Without the porous 
plug, there would be no significant fountain pressure produced, because the 
channel diameter is too large. In this case, the pressure within the heat pipe 
would be given by point b', which is within the vapor region of the phase diagram 
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and would result in cavitation and vapor lock of the loop. 


A third possibility is to have the porous plug present at the terminations of the 
loop, but to have larger pores so that there is some fountain pressure achieved but 
it is not as large as predicted by London’s equation. The advantage of larger pores 
is in the plugs’ improved ability to transport heat from the loop to the storage 
dewar. The important characteristic of the intermediate pore size model is to 
ensure that the pressure rise, due to the fountain effect, be larger than the slope 
of the vapor pressure curve, 


f dP 
^dT 


> 


^dP N 
dT 


>vp 


(7.14) 


Combining eqs. (10) and (12), we get that 




dP 

dT 


A 




1 


ps 



(7.15) 


Also, using the Clausius-Clapeyron equation and assuming ideal gas behavior we 


get 
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_ X _ Xp 

v dTjsvp _ Tv g RT 2 


(7.16) 


where X is the latent heat of helium and v g is the specific volume of the vapor. 
The above two equations can be combined to get the conditions in heat flux and 
pore diameter in the plug 

qd < 



( Y 

g 

ps RT 2 _ 1 

f 

Xp 1 


1/2 


(7.17) 


Inserting values for ASTROMAG at 1.8 K, q = 500 W/m 2 requires that d < 1.3 mm 
which is quite large for porous media. On the other hand, during charging, the 
heat flux will increase by a factor of 25 which would therefore require d < 50 pm. 
Listed in Table 5 are the maximum diameters for a variety of operating conditions 
and temperature as calculated from eq. (7.17). To ensure sufficient heat removal 
and safe non-cavitating operation of the heat pipe, the plug should have a pore 
size around 10 pm. It also would be advantageous to fabricate the plug of a high 
conductivity metal, such as copper, to enhance the heat exchange between the 
helium within the loop and the bath. 

The use of a He II heat pipe for cooling of ASTROMAG would seem to be justified 
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based on the above arguments. As envisioned, the system would consist of a 
cooling loop containing static He II terminated by porous plugs immersed in the 
storage dewar. Steady state and charging heat leaks could easily be removed by 
this system. However, the question remains as to the behavior of the system 
during extraordinary operation such as initial cooldown or magnet quench. These 
modes of operation will require active control. 

First, consider the behavior of the He II heat pipe during magnet quench. One of 
the most important issues here is to minimize the loss of helium inventory. Since 
the He II heat pipe contains relatively little helium of its own (approximately 10 
liters), this helium would need to vent due to rapid vaporization. Relief values 
(RV) are required for this purpose (see Fig. 19). Subsequent to the quench, the 
magnet is expected to be at relatively high temperatures (100 - 120 K). Cooldown 
and resupply would therefore be accomplished by opening the vent valve to space. 
This process would allow helium to circulate through the loop. Initially, the loop 
would be filled with low pressure vapor, thus cooling would be rather slow. In 
fact, the porous plug might act as a phase separator in this situation preventing 
the liquid from entering the loop. To ensure rapid cooldown, the porous plug 
would need to be equipped with a heater to keep the downstream temperature 
above that of the helium within the storage dewar, thus activating the fountain 
effect. Helium would be circulated through the loop until the magnet is returned 
to operating temperature. Initially, the helium would vent, but subsequently, 
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when He II is being circulated. Opening of v 2 would allow recovery of the helium 
to the storage dewar. A similar procedure would apply when cooling the magnet 
from initial ambient temperature. 


Using the porous plug as a fountain pump should be straightforward since the 
porous plug is intended to establish a significant fountain pressure during steady 
state operation. However, because of the need to extract heat through the plug, it 
would not be expected to operate as an ideal fountain pump. Even for an ideal 
pump, the mass flow rate is reduced below that given by eq. (7.3) as the result of 

i q 

hydrodynamic pressure loss. In general, the mass flow rate is given by 


m = 


$ 

6.6 * s in T in 

P 


( 7 . 18 ) 


For high mass flow rates such as are required for SHOOT or may be required for 
ASTROMAG during extraordinary operation do result in sizable pressure drops 
and therefore reduced mass flow rates. Also, because of the heat leak due to 
normal fluid flow through the porous plug, the heat required to achieve a given m 
would be larger than predicted by eq. (7.18). 
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7.6 Conclusion 

A He II heat pipe based on the principals described in this report should work 
effectively for heat removal from the magnets for ASTROMAG. The primary heat 
exchange mechanism would be that of internal convection. As envisioned, the 
heat pipe would consist of a tube wrapped around the magnet having either end 
ter mi nate within the helium storage dewar. Each end of the tube has a porous 
plug which serves to prevent cavitation while at the same time allowing heat 
transfer. Based on preliminary calculations, a plug having the same cross-section 
as the tube but consisting of multibore tubes of diameter 10 pm should perform 
adequately. Further analysis and experimentation would be required to verify 


these findings. 
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Table 5. Maximum pore diameter to prevent cavitation. Ap = 0.4 cm , Q = 10 mW and 250 mW, a = 0.5. 
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8. HOFMANN’S LOOP SYSTEM 


8.1 Introduction 

The central idea of Hofmann’s Loop System (HLS) is to use object’s anyway 
dissipated heat for driving a Fountain Effect Pump (FEP) and sustain an external 
convection He II cooling system. Furthermore, one might expect it to be a self- 
sustained cooling loop in the sense of feedback system which means: a transient 
increase in heat load will drive stronger the FEP, increase the flowrate and its 
cooling capacity. Decrease in dissipation will affect in the opposite direction. It 
was suggested and analyzed theoretically in source (14) and checked through a lab 
model as summarized in sources (15), (16) and (17). 

However, some difficulties arise while trying to apply that idea for cooling the 
ASTRQMAG. In the previous chapter (7) a comparison was drawn between 
internal and external convection. Since the HLS is one of the leading design 
concepts the following chapter is a more detailed attempt to derive its main 
parameters to meet the ASTROMAG’s cooling and performance specifications (2,3) 
for demonstrating its shortages. 

8.2 Non-Dissipated Heat Load 

While the dissipated energy of the magnet is used to drive the FEP, part of it "is 
saved" from being dissipated as mechanical work. It is the work for circulating 
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the He II fluid and does not account as boil-off. 

However, an idealized analysis based on perfect heat exchangers (14) shows that 
the maximum energy porition saved and converted to mechanical work does not 
exceed 10%. That is consistent with a result obtained in (13) saying that the 
maximum efficiency of a FEP is about 15%. 

Practically, with non-ideal heat exchangers one should expect about 5% of the 
dissipation saved as mechanical work. 

8.3 Flow Velocity and Mass Flowrate 

Matching the ASTROMAG steady state requirements was demonstrated in 
previous chapters through conduction of thermal straps. The equivalent thermal 
conduction of He II (or its internal convection) is higher by about two orders of 
magnitude than that of any metal. So, just the internal convection is adequate for 
the steady state cooling. Therefore, the criteria for minimal HLS velocity is 
suggested: the heat transport through convection should be about twice the 
conductive heat transport. Put another way the total heat transport is required to 
be about 3 times the conductive (internal convection) heat transport. 

Van Sciver (11) showed for small AT that external and internal heat transport 
relation is ruled by a non-dimensional parameter K, defined as: 
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K - p . c p . V . (fL) 1 ' 3 fT m - T b ) 2 ' 3 


where: v - velocity 

f - temperature dependent He II heat conductivity function (source) 
L - channel length 

T m - T b - temperature difference across the channel 

To match the above criteria for minimal HLS velocity, the non-dimensional 
parameter K has to reach the value of 5. 

T = 1.4 -r 1.5 K 
p = 0.14 G • cm' 3 
c p = 1 J • g 1 • K' 1 
f(T) = 1/300 cm 5 • K • W' 3 
Tm ’ T b = 0.3 K 
K = 5 


According to source (3) the length of the channel might be 20 to 40 m and its 
diameter 6 to 8 mm. 


For: L = 20 m and D = 8 mm. 
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if D = 0.6 cm then m = 2.14 g/s 
For L = 40 m and D = 0.6 cm 
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v = 33.6 


cm 

s 


m = 1.7 1 

S 


For L = 30 m and D = 0.6 cm 


v = 37.0 


cm 

s 


m = 1.86 JL 

s 


8.4 FEP’s Driving Heat Power 

The interaction of a FEP and a transfer line was formulated and solved in source 
(13). Following that procedure some heat powers and pressure drops are derived 
for L = 20, 30, 40 m and D = 6, 8 mm, T = 1.4, 1.5 K. For small velocities f is 
approximated as 0.004 to simplify the calculations. 
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8.4.1 L = 20 m 


T = 1.4 K 


D = 6 mm 


f= 0.004 


Q (W) 

0.05 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

1 

m (g/s) 

0.27 

0.53 

1.02 

1.45 

1.82 

2.14 

2.43 


AP (Pa) 

0.017 

0.065 

0.24 

0.48 

0.76 

1.05 

1.35 

ES 


AP m ax = 70 Pa 

Qmax = 77.22 W 

“max = 229 
(f= 0.00231) 


8.4.2 L = 30 m 


T = 1.4 K 


D = 6 mm 
f= 0.004 


Q (W) 

0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

m (g/s) 

0.27 

0.52 

0.99 

1.38 

2.00 

2.48 

3.02 

AP (Pa) 

0.025 

0.095 

0.34 

0.66 

1.39 

2.12 

3.16 


4 P max = ™ Pa 
Qmax - 59 4 W 

"W = 181 e/* 
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8.4.3 L = 40 m T = 1.4 K 

D = 6 mm 
f= 0.004 


Q (W) 

0.05 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

m (g/s) 

0.27 

0.54 

1.05 

1.53 

1.95 

2.35 

2.70 

3.00 

AP (Pa) 

0.008 

0.03 

0.12 

0.25 

0.42 



|d 


AP max = 70 Pa 
Qmax = H2.9 W 

"W = 335 ef* 

(f= 0.00231) 


8.4.4 L = 20 m T = 1.5 K 

D = 6 mm 
f= 0.004 


Q (W) 

0.05 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

m (g/s) 

0.17 

0.34 

0.67 

0.99 

1.29 



2.07 

2.3 

2.51 

2.70 

AP (Pa) 

0.007 

0.026 

0.10 

0.22 

0.38 



0.99 

1.21 

1.45 

1.68 


APmax = 68 p a 
Qmax = 76 4 W 
mmax = 22.6 g/s 
(f= 0.0023) 
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8.4.5 L = 30 m 


D = 6 mm 


T = 1.5 K 


Q (W) 

0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

m (g/s) 

0.17 

0.34 

0.66 

0.97 

1.51 

1.97 

2.54 

AP (Pa) 

0.01 

0.04 

0.15 

0.33 

0.80 

1.35 

2.22 


APmax = 68 Pa 
Qmax = 61-3 W 

I * 1 max = 18-1 


8.4.6 Conclusion 

Let us refer to the case of L = 30 m and D = 6 mm as a representing one. For 
T = 1.4 K a heat power of about 0.45 W should be applied for reaching the 
1.86 g/S flowrate. At T = 1.5 K it increases to about 0.65 W. Those required heat 
power supplies have to be compared to the (assumable) specified magnet 
dissipation: 0.01 W per coil. Comes out the conclusion that there is not enough 
driving heat power for supporting a significant external convection, and it becomes 
an internal convection case. 


If we would consider an external supply of about 0.5 W to support that flowrate 
than we increase the total heat load and He II consumption by about one and a 
half order of magnitude which is not acceptable either. 
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8.5 Stability 

The self-sustained FEP operation is a closed loop feedback system. As such it 
should be examined about the conditions of stability. 

Hofmann describes^ 17 ^ the first unsteady heat load experiment as undergoing 
strong oscillations. The pressure oscillated between zero and its maximum value 
by a time constant of about 150 especially in the low values of FEP driving heat 
powers. Hofmann reports that in order to ensure stable behavior an additional 
constant heat supply of about 0.5 W is required for complete suppression of the 
oscillations. That additional heat power seems to be incomparable to the 0.02 W 
magnet steady state dissipation. 

A block diagram in Fig. 21 is suggested for system engineering analysis of the 
HLS. On its basis a stability condition might be derived for non-oscillating 


behavior. 
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Fig. 21: Block diagram linearized representation of Hofmann’s Cooling Loop 

for system stability analysis. 
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